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Аннотация 

Введение. Функционально-градиентные материалы находят большое применение, т.к. неоднородность свойств 
позволяет управлять прочностью и жесткостью контрукций. Этим вызван большой интерес к данной теме в 
мировой научной литературе. Построение решения таких задач существенно зависит от типа граничных условий. 
В настоящей работе рассматривается равновесие тонкостенного кругового цилиндра, механические свойства 
которого заменяются вдоль радиуса. На цилиндрических поверхностях заданы однородные граничные условия, 
которые до этого не рассматривались, воздействие оказывается на торцах. Математическая постановка задачи 
осуществляется в линейной теории упругости в рамках осесимметричной деформации. В работе построены 
выражения для компонент напряженно-деформированного состояния цилиндра, в которых некоторые 
коэффициенты находятся из решения полученной системы линейных алгебраических уравнений. 

Материалы и методы. Материал цилиндра является линейно упругим, модуль упругости которого линейно 
зависит от радиальной координаты. Основным методом исследования является ассимптотический метод, в 
котором в качестве малого параметра выступает половина логарифма отношения внешнего и внутреннего 
радиусов. Для построения характеристик напряженно-деформированного состояния цилиндра применены 
итерационные процессы. 

Результаты исследования. Для линейно-упругого функционально-градиентного полого тонкостенного 
цилиндра получены однородные решения краевой задачи. Анализ этих решений позволяет раскрыть характер 
напряженно-деформированного состояния в стенке цилиндра. С этой целью проведен асимптотический анализ 
решений, получены соотношения для перемещений и напряжений. Установлено, что эти решения соответствуют 
пограничному слою, при этом их первые члены определяют краевой эффект Сен-Венана, аналогичный теории плит. 
Обсуждение и заключение. Построенное с помощью асимптотического разложения аналитическое решение 
задачи о равновесии неоднородного по радиусу тонкостенного цилиндра может быть использовано для 
численного решения конкретной задачи. Для этого нужно решить полученные системы линейных 
алгебраических уравнений и определить соответствующие коэффициенты. Полученные асимптотические 
представления позволяют анализировать трехмерное напряженно-деформированное состояние. Выбор 
количества членов разложения позволяет рассчитать перемещения и напряжения с заданной степенью точности. 
Этот анализ может быть полезен при оценке адекватности прикладных методов расчета, применяемых в 
инженерной практике. 


Ключевые слова: линейная теория упругости, функционально-градиентный материал, тонкостенный полый 
цилиндр, однородные решения, пограничный слой, вариационный принцип 
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Введение. Функционально-градиентные материалы находят широкое применение в различных конструкциях. 
Благодаря зависимости механических свойств от координат можно управлять напряженно-деформированным 
состоянием (НДС) деталей. Примером использования такой неоднородности может служить цилиндр, 
механические свойства которого зависят от радиуса. При этом интерес может представлять цилиндр как 
отдельная конструкция, так и являющийся подтелом составного тела, например, соединяющий две среды с 
сильно различающимися свойствами. При расчете НДС тонкостенного цилинда используются некоторые 
прикладные теории. При этом важно знать оценку их адекватности. Особенно в случае неоднородных свойств 
она может быть осуществлена с помощью компьютерного моделирования или асимптотического анализа на 
основе трехмерной постановки. Последним определяется актуальность настоящего исследования. 
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Изучению НДС полых цилиндрических тел в рамках линейной теории упругости посвящен ряд исследований. 
В [1, 2] на основе метода сплайн-коллокации и метода конечных элементов изучено механическое поведение 
радиально-неоднородного цилиндра в трехмерной постановке. В [3] изучено НДС цилиндра, свойства которого 
зависят от радиуса, нагруженного равномерным внутреннем давлением. В [4] проведено аналитическое 
исследование для функционально-градиентного пьезоэлектрического цилиндра. В [5] построено точное решение 
для радиально-неоднородного полого цилиндра с экспоненциальным модулем Юнга, постоянным 
коэффициентом Пуассона и степенным модулем Юнга. В [6,7] с помощью метода прямого интегрирования 
получено аналитическое решение осесимметричной задачи термоупругости для сплошного цилиндра, когда 
коэффициент линейного теплового расширения является произвольной функцией от радиуса. В [8] разработана 
общая асимптотическая теория трансверсально-изотропного однородного полого цилиндра. Для трансверсально- 
изотропного однородного цилиндра получены новые группы решений. Приведено сравнение построенных 
решений с решениями, построенными с помощью прикладных методов расчета. В [9, 10] изучены некоторые 
краевые задачи теории упругости для функционально-градиентного изотропного и трансверсально-изотропного 
(плоскость изотропии перпендикулярна оси) цилиндра, в случае, когда модули упругости являются 
произвольными непрерывными функциями от радиуса цилиндра. В [11] проведен анализ задачи изгибной 
деформации для радиально-неоднородного цилиндра. Анализ этих работ показывает, что не для всех типов 
граничных условий на цилиндрических поверхностях имееются асимптотические представления решений. 

В настоящей работе на основе асимптотического анализа трехмерных уравнений теории упругости изучаются 
особенности НДС тонкостенного цилиндра, свойства которого линейно меняются вдоль радиуса. При этом 
внутренняя граница закреплена в осевом направлении и свободна в радиальном направлении. 

Достижение этой цели было получено с помощью нескольких шагов: асимптотическое интегрирование 
дифференциальных уравнений и построение однородных решений; вывод формул для компонент вектора 
перемещений и тензора напряжений; учет краевых условий на торцевых поверхностях. 

Материалы и методы. Рассматривается радиально-неоднороднополый тонкостенный цилиндр 


Г= { Е в фе [0; 2*|, = [-4;% ] в цилиндрической системе координат г,ф, < с началом на его оси. Задача 


о его равновесии, в случае закрепления цилиндрических поверхностей вдоль оси и нулевых нормальных 
напряжений, решается в осесимметричной постановке при действии напряжений на его торцах. 
Краевая задача состоит из уравнений равновесия [8]: 


ых ха +— 2 =0, (1) 
ди 05 г 
05. + 05. + 0, = 0, (2) 
дг 05 г 
где о„,о„,О,б„— компоненты тензора напряжений. 
Определяющие соотношения [8]: 
а. = 6+4) 7% [9 (3) 
дг г 
с. - (26+%) 9% +29 = (4) 
ы 05 д г 
и ди. ди 
ен =(2“+^)— +7 —+— |, р 
т ) г 0 0 © 
с. = с“ + = (6) 
. 0% бк 
Здесь и, =и,(г,5), и, =и, (г, <) — компоненты вектора перемещений. 
Параметры Ламе изменяются по линейному закону вдоль радиуса: 
С(г)=С.т, А(г) =, (7) 
где сС.,Л. — константы. 
После подстановки (3)—(7) в уравнения (1), (2) безразмерная система уравнений принимает вид: 
0?и ди, д?и. ди. 0?и, 
(26, +) +8—2 |+=(б, +2, ее — = +А.Е?е® — + Су ?е2® —_^ — 2СуЕ?и, = 0, (8) 
др? д дрдё 06 08? 
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ди, би. б?и. ди д?и 
С, +—* |+(26, +, )=?| е?® че —^ |+(б,+^, вет —2=0. (9) 
др? др ОЕ? 08 дрдё 
Здесь: 
1 т м 1 Г, ‚. 
р=-Шш| —|, &=— — новые безразмерные координаты; &=-—Ш|-— | — малый параметр в случае 
5 п ”® 2 п 
| А. С. 
тонкостенности; п =./, ре [181], 6 =[-[;1, 1=—°; и, = и, ‚ = а Аве ло ‚ б= ло ; С 
® п п С, С, 


некоторый параметр, имеющий размерность напряжения. 
Рассмотрим задачу, в которой на боковых поверхностях цилиндра заданы однородные граничные условия: 


и: =0, (10) 
би =0. (11) 
К торцам цилиндра приложены напряжения: 
бы =» (В), (12) 
би =, (р), (13) 
(+=62) 
[ [} 


УГ 72 ыы 55 


с = с’ 6 = с’ б= С безразмерные напряжения. 
1 1 1 


Компоненты вектора напряжений 1, (р),6,, (р), (5=1;2) удовлетворяют условиям равновесия. 


Для построения однородных решениий компоненты вектора смещений будем искать в виде: 


и, (р;5) =и(р)е“е, и (р;&) = и(р)е*. (14) 
Подставляя представления (14) в систему (8)—(11), получим: 
(26, +^,) (и" (р) + ги’ (р))+ ве ((б,+7, м (р)+.=и(р))+ =2С, (а?е?> —2)и (р) = 0, (15) 
5. (и"' (р) + вм" (р)) +(2С, +^., )=? (аети (р) + че? (р) +8(С,+^. )аети’ (р) = 0, (16) 
У. =0, (17) 


= 0. (18) 


р=-+1 


[(26, + м (р)-+=^. (и(р)+ чету (р)) | 


Исследуем краевые задачи (15)-(18) при &—0. Для решения (15)—(18) при &-—0 воспользуемся 
асимптотическом методом [9-13]. 
Ненулевые решения (15)—(18) соответствуют третьему итерационному процессу, компоненты вектора 
смещений будем искать в виде разложений по малому параметру: 
и(3) (р) = =(и (р) + 5И.1 (р) я .), 


из (р) ==(и (р)-+ =, (р)+---), (19) 
а = =! (В +В, +---). 


После подстановки разложений (19) в уравнения (15)—(18) для членов первого порядка имеем: 


(26, +2 (р)+Вь (6, +2) (в) бы (р) =0, со) 
Сом (р) + В, (С, + изо (рР)+ (26, +0 Ви (р) =0, (20 
У (Ри. =0, (22) 


= 0. (23) 


р=+1 


(26+ )мь (р) + Вим (р) 


Следуя [13], спектральная задача (20)—(23) соответствует потенциальному решению для плиты. 
Таким образом решения представляются в виде: 


а) ибо (р; &) = =>. т, (-26. В, эпВих $ (Вор) = 0(=)) й оч (Ви +В, +.- .) :] ь (24) 
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> : 1 
м" (р; 5) = => (2.В, зшВь, со$ (Вир) + 0(=)) * вар (ры +аВ,, +.- .) :) . (25) 
Здесь В.„ является решением уравнения: 
соз Вх = 0. (26) 
Напряжения, соответствующие решениям (24) и (25), имеют вид: 
ов = т (- —4С2Вз, зтВих с0$(Вир) + О (Е (одуееме[ Вох + ЕВь +8) , (27) 
0. = => Т, (4 ?В, п Во т (Вир) + О(= х еек Вик +=Вь, +-- 3) ь (28) 
& 
е 1 
бе 9 = =>, (4 СоВа 91п Вох соь(Выр)-+0())екр] (ры +В, +=) (29) 
5 = 0(:). (30) 
[9] 1 
Ь) из (р;Е) = =>. :(2С.В:, созВ., соз (Вир) +О(= =))х ор +В, +.. 3} (31) 
Э Ь : 1 
иг” (р; 5) — => Е (26,88, созВь, п (Вир) ы: 0(=)) 2 ее: (Во +=В,, +-- .) :] . (32) 
Здесь Во, является решением уравнения: 
зшВь; = 0. (33) 
Напряжения, соответствующие решениям (31) и (32), имеют вид: 
[2] . т 
0%? =>. Е (-403В% созВо (Вир) + ое +аВ, +. 3} (34) 
со 1 
6%? = =>Е ‚ (4С2В, созВо, соз (Вир) +О (= ево + В, +] (35) 
Ьы : 1 
од = 52 (4010 воз, м (Вьр)+ (6) вы ару += в6) 
68? =0(=). (37) 


Общее решение (15)—(18) будет суперпозицией решений (24), (25), (31), (32): 
и, (р;5) = =>. Т, (-26.В&, зпВо, т (Вир) + 0(=))х 


С] 


ке] (ры +В» ++ =>. Е :(2С,В&, созВ; соз(Выр) + обе] (р + &В,, + =) (38) 


и, (р;5) — =>, (26,В, зп Ви соз(Вир) + 0(=))х 


кем (Вы + ёВ, + ...) в (26,В3, созВ, чт (Вор) +О(=)х окр + В, + ...) :] (39) 


Решения (24), (25), (31), (32) имеют характер пограничного слоя. При удалении от торцов решения (24), (25), 
(31), (32) экспоненциально убывают. 
Чтобы определить константы Т., ЕЁ, воспользуемся вариационным принципом Лагранжа. Вариационный 


7 


принцип принимает вид [8]: 


21 
> ] [(с„. 1, ) би, + (5 —ь,, би. | е2 р = 0. (40) 
— = 
Подставляя (24—36) в (40), имеем: 
Х Мать = Ро, 41 
К=1 
Зо, в) о 
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Здесь: 


М „ =16б0 ВВ (Во, —Ва» ) зп Вих Во, чп (Во, — В» )х =з- А беты › (при а к) 
‘ = = 


, 28, 1 ов, 
М ; =16С3В8, п? Ву, она (при =К) 


| ы 
Ро = 2СВ5, яп Во; | (> (р)соз (Выр)- и бррр уар-ече е } 


+1 (6. (р)есз (Вир) и. буть очен 


г ог +Во; )1 + Во; |1 
О; =166. В Ви (Во В, ) созВь, ; созВ,, чл (Во; — В» ) эВ ыы ‚ (при #* 1) 
, - | 
О; =16С5Вь, со? Во, |= т ны (при 1=] ) 
| = & 


(ею (Выр) + бумер }ар-ено м, 
7. 


Постоянные Т,,Ё, (р=1,2,...) находятся из систем линейных алгебраических уравнений (41), (42), 


Ру = 2С.Вь; созВь, | 


Во 


7 
= 


+ (6. (р)соз(Вор) + (р)з (Выр))ар . 5 


Т, =То +ЕТи+.-*, 
Е =ЕВ +ЕЁ, +... 


аналогичные которым изучены в [13]. 

Результаты исследований. В работе в осесимметричной постановке рассмотрено решение задачи линейной 
теории упругости для функционально-градиентного полого тонкостенного цилиндра, свойства которого 
изменяются по толщине по линейному закону. На боковых поверхностях цилиндра заданы однородные 
перекрестные граничные условия, на торцах задан вектор напряжений. Построенные однородные решения 
удовлетворяют граничным условиям на цилиндрических поверхностях. Для их построения использован 
асимптотический подход на основе раздожения по малому параметру, характеризующему относительную 
толщину цилиндра. Для учета неоднородных граничных условий на торцах получены системы линейных 
алгебраических уравнений, аналогичные изученым в литературе. Показано, что построенные решения НДС 
имеют погранслойный характер, который соответствует краевому эффекту, аналогичному теории неоднородных 
плит, носящему имя Сен-Венана. 

Обсуждение и заключение. Обычно при изучении НДС тонкостенных конструкций строятся прикладные 
методы расчета, снижающие размерность задачи. В этой связи актуальной является задача определения 
диапазона геометрических и механических параметров, в которых эти методы дают приемлемую точность. 
Построенные в работе на основе асимптотического анализа решения трехмерных уравнений позволяют оценить 
адекватность таких прикладных теорий с заданным наперед порогом точности. Кроме того, эти решения могут 
найти применение при оценке численных решений задач для конструкций с функционально-градиентными 
материалами. 
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